Decay rates for topological memories encoded with Majorana fermions 
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Recently there have been numerous proposals to create Majorana zero modes in solid state het- 
erojunctions, superconducting wires and optical lattices. Putatively the information stored in qubits 
constructed from these modes is protected from various forms of decoherence. Here we present a 
generic method to study the effect of external perturbations on these modes. We focus on the case 
where there are no interactions between different Majorana modes either directly or through inter- 
mediary fermions. To quantify the rate of loss of the information stored in the Majorana modes we 
study the two-time correlators for qubits built from them. We analyze a generic gapped fermionic 
environment (bath) interacting via tunneling with different components of the qubit (different Ma- 
jorana modes). We present examples with both static and dynamic perturbations (noise), and using 
our formalism we derive a rate of information loss, for Majorana memories, that depends on the 
spectral density of both the noise and the fermionic bath. 
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I. INTRODUCTION 

Topological quantum computation requires the exis- 
tence of topologically ordered states whose low energy 
excitations follow non-Abelian statistics. The subspace 
of states corresponding to a fixed number of quasipar- 
ticles is degenerate, to an exponential precision, in the 
separation between quasiparticles, and an exchange of 
the positions of these anyonic excitations, also known as 
braiding, leads to a unitary transformation within this 
low energy subspace. These unitary operations are insen- 
sitive to the exact path used to perforin the braiding op- 
eration and in many cases, for an appropriate encoding, 
braiding operations correspond to "standard"' one- and 
two-qubit gates within the low energy subspace. These 
operations can be used as building blocks for fault toler- 
ant quantum computation. 

There are many candidate systems for experimental re- 
alizations of topological phases of matter with these prop- 
erties. There is preliminary evidence that the v = hjl 
fractional quantum Hall state may have non-Abelian 
excitations^"— . Spin-triplet p^^ -I- i'Py pairing superfluid- 
ity occurs in the A-phase of "^Heii^ and in strontium 
ruthenates^"— , in which half quantum vortices would 
be non-Abeliani^ii^. There are also proposals to re- 
alize chiral p-wave superconductors in ultra-cold atom 
systemsi^"— . Furthermore there have been many ad- 
vances towards producing topological states of matter in 
layered heterojunction systemsii"— . 

Virtually all current experimentally viable proposals 
for platforms for topological quantum computation only 
support Ising type anyons which are carried by Majorana 
fermion modes. Colloquially speaking these fermions are 
half of a regular fermion. More precisely they are self- 
adjoint operators 7^ which can be written as a sum of an 
annihilation and creation operator for one fermion mode 
and which satisfy the algebra: 



{7,, 7j} = 2%, 7/ =7i 



(1) 



into a regular fermion mode c and its adjoint c^ via c = 
\ (71 + «72) and c^ = \ (71 - ^li)- 

The topological qubit is made up of four spin polarized 
MBSs 71, 72, 73 and 74^. These can be combined into 
two sets of creation and annihilation operators: 

(2) 
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For the logical basis it is convenient to work in the even 
fermion parity subspace. The qubit basis can be chosen 
to be |-fi,) EE |00) and |-l) eee |11) where the O's and 
I's refer to the occupation numbers relative to the com- 
plex fermion operators in Eq. ([2]). Because of fermion 
parity conservation, any operation that does not entan- 
gle the states with the environment cannot mix even and 
odd fermion parity states for the qubits. As such, all 
gates acting on the topological qubit should not take the 
system out of the logical subspace. Furthermore all the 
operators of the single spin Clifford group may be pro- 
duced by braiding the four vortices of our qubit leading 
to potentially topologically protected gates^^. In partic- 
ular the various single qubit operations in our logic basis 
may be conveniently written in terms of the Majorana 
operators. For future use we note that in this encoding 
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Any two Majorana fermion operators can be combined 



Here all the sigma matrices are with respect to the 
logic basis |-I-l) and |— l). We will primarily be in- 
terested in correlators of the form (cr^ (0) cr^ (T)) = 
- (71 (0) 72 (0) 71 (T) 72 (T)). We will proceed to calcu- 
late these below. 

The Majorana operators are zero modes of some mean 
field Hamiltonian [i^MF, 7i] = so it can be argued that 
these modes are protected from decoherence as the mean 
field Hamiltonian when restricted to the subspace gen- 
erated by these modes is zero. One of the open tasks 
of topological quantum computation is associated with 
understanding the extent of this protection. This is the 
subject of this paper. 



II. SUMMARY OF MAIN IDEAS 

In this section we outline the setup of the rest of the 
paper. We present the relevant Hamiltonian and discuss 
its basic properties. We describe the type of qubit we will 
focus on in the text, a localized Majorana mode, and give 
an overview of some other encodings we shall not consider 
in this paper. We describe the kinds of calculations of 
memory coherence we are going to do in this paper. We 
also give a Section by Section outline. 

We begin our discussion with relevant Hamiltonians. 
The Majorana fermions interact with the external envi- 
ronment via tunneling type Hamiltonians. On symmetry 
grounds, for a single Majorana mode, any such interac- 
tion may be written as: 

Hint ^1 f d'^r [uoir) $^ (r^ *^ (f) - uS(f) * (r) $ (t^] . 

, , (4) 

Here uq (r) is the localized mode function associated with 
the Majorana bound state, $(r) is any local bosonic field, 
which in the simplest case is a tunneling amplitude (com- 
plex number) and ^(f^ is a regular (complex) fermion 
field. In this paper we will analyze multiple Majorana 
fermions coupled to different types of environments via 
Hamiltonians of the form given in Eq. Q . Furthermore 
the fermions in the bath will always be assumed to be 
gapped, for example, electrons in an insulating or super- 
conducting material (environments composed of gapless 
fermions, instead, would obviously lead to decoherence) . 
There are many examples of microscopic situations 
where Hamiltonians of the form given in Eq. ^ arise, 
one is as follows. If one writes the mode expansion of 
the electron creation and annihilation operators in the 
(superconducting) system of interest, one finds that: 
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Here qe stands for the eigenoperators of the BdG equa- 
tions, with non-zero energies, while ue and ve are the 
components of the corresponding eigenmode of the BdG 
equations. 7 is the Majorana fermion corresponding 
to the zero energy mode. Now consider an insulat- 
ing substrate below a system which may be described 
by Eq. ([5|) above. A concrete example is given by 
the bulk of a topological insulator in tunneling contact 
with a superconductor as shown in Ref. l33l . For a 
static Hamiltonian the bulk and surface states are or- 
thogonalized, but dynamical effects such as phonons or 
two-level defect systems can alter the original Hamilto- 
nian and turn on a hybridization. This perturbation 
takes the form of a tunneling between the electrons: 
Hint = J d'^r ^ (r) *^ (f) * (r) + h. c, where $ (f) con- 
trols the amplitude of fluctuations of the tunneling cou- 
pling. $ (7^ can be due to phonons, two-level systems, or 



even classical sources of noise. The electrons ^ (r) come 
from the insulating (gapped) system, which comprise the 
fermionic component of our bath. This illustrates one 
of the many ways to arrive at Hamiltonians of the form 
Eq. Q. 

The coupling Hamiltonian that is derived in the para- 
graph above is local. The terms in Eq. (j4]) are local and 
couple to only one Majorana mode, with no long dis- 
tance coupling between the modes of any form. In this 
paper we shall focus on sets of baths that couple to each 
Majorana individually. We would like to stress now and 
henceforth that even by coupling to individual modes, 
one at a time (with no cross mode coupling), the bath 
can be very damaging, in many cases leading to zero co- 
herence for long times. 

Below, we look at decoherence by analyz- 
ing qubit correlations such as (cr^ (0) a^ (T)) = 
— (71 (0) 72 (0)71 (T)72 (T)), which, as we show in 
this paper, factorizes when the baths that couple to each 
Majorana are uncorrelated with one another: 

{a^ (0) a- (T)) = (71 (0) 71 (T)) x (72 (0) 72 (T)) . (6) 

Thus, even though the qubit is defined non-locally us- 
ing spatially separated Majorana fermions, below we will 
show that the decay of the memory is controlled by the 
product of the two-time correlations of the separate Ma- 
jorana modes. It then suffices to understand the effect of 
the bath on each Majorana fermion separately. 

At this point its worthwhile to stress that the qubit 
encoding given above is not unique. A particularly 
interesting example of a different encoding, given by 
Akhmero\s^, is a fermion parity protected encoding. 
There, the qubit is made from fermion parity preserv- 
ing operators: 



7 = 7 
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that commute with both the tunneling Hamiltonian and 
the Hamiltonian for the environment. Here the Cj are the 
operators in the mode expansion of the fermionic ^ (r) 
field in the bath (i here labels the mode, which can be mo- 
mentum, for example). For a finite system, such as mid 
gap Carroll Matricon deGennes states in vortex cores, 
this compound qubit is very efficient. However we stress 
that, in the presence of a bath (say made by continuum 
states) , the construction of an operator that is protected 
because of parity conservation requires a product of in- 
finitely many operators: which is not practical or easily 
experimentally measurable. One could also truncate the 
product so as to account for a system, and the terms 
omitted are those assigned to the bath, as depicted in 
Fig. [T] In this case, however, because the operator lacks 
degrees of freedom assigned to the bath, parity can leak 
to the environment decohering the qubit. As such we 
will ignore all "compound" encodings for the rest of the 
paper. 

Finally, we would also like to mention that the above 
scheme, with simple, non-compound, Majorana encod- 




Figure 1: Depiction of the separation between system and 
batii degrees of freedom. For infinite batiis, one cannot con- 
struct a local operator of the form Eq. ([7|, one that is a prod- 
uct of a finite number of terms. If the product is truncated, 
parity leaks into the bath. 



ing, generalizes to multiple qubits. One possible encod- 
ing (though not the most economical) is to use four vor- 
tices and as such four Majorana modes per qubit. For 
this and any other encoding all possible correlators for 
the quantum memory may be expressed as expectation 
values of various products of Majorana operators^. All 
quantum coherences for our qubits may then be com- 
puted by studying Majorana mode correlators which we 
study below. 

In carrying out this program, we will analyze two dis- 
tinct types of environments: the first is when couplings 
$(r) change suddenly but remain static thereafter, and 
the second when the environment changes dynamically. 
We show that that in the static environment case the tun- 
neling Hamiltonian merely leads to a finite depletion of 
the Majorana two-time correlations. In this case, much 
of the information stored in these modes survives for ar- 
bitrarily long times. 

More generally, for dynamic environments, we obtain 
an expression for the rate of loss of information stored 
in the Majorana operators that depends on the spectral 
density of the noise and of the fermionic bath. We present 
several examples of noise that can be studied essentially 
exactly, for instance classical telegraphic noise, as well as 
both classical and quantum Gaussian fluctuations. 

The results in the paper are presented as follows: 

• In Section IIIII we present general considerations 
involving the coherence properties of Majorana 
modes. We show that under reasonably generic 
initial conditions the coherence of the Majorana 
modes does not depend on their initial states. Fur- 
thermore we show that the two time correlation 
functions, coherences, factorize as a product over 
coherences for individual Majorana modes, that 
make up the quantum memory, interacting with 
their individual environments. As such we may re- 
duce the problem of the coherence of the quantum 



memory to the problem of the coherence of one Ma- 
jorana mode in tunneling contact with a (gapped) 
fermionic reservoir. 

• In Section ITVl we take a first step towards a calcu- 
lation of the coherence of a single Majorana mode. 
We begin by describing the Keldysh technique rel- 
evant to Majorana modes. We present combinato- 
rial tricks that make is possible to efficiently con- 
vert Keldysh computations using a mixture of Ma- 
jorana and regular fermionic modes into a more fa- 
miliar computation which uses only regular fermion 
modes. We then present an example where, for 
simplicity, we treat the fermions in the bath as free 
(non-relaxing approximation). We also present a 
general formula for the coherence of a Majorana 
qubit that is used several times in the remaining 
analysis. 

• In Section |V] we present several related classical 
models for the fiuctuations of the bath. We solve 
these models essentially exactly, by mapping the 
problem of the coherence of a single Majorana mode 
to the problem of a particle undergoing classical 
diffusion. We use this technique to study classical 
fluctuations of the tunneling amplitudes and energy 
levels of the reservoir (we primarily focus on Gaus- 
sian fiuctuations). In all cases we find decoherence 
with a rate that depends on the spectral density 
of the fluctuations in the reservoir. In many cases 
the decoherence due to an individual fermion mode 
has a power law time dependence but it will turn 
out that a bath made of many weakly interacting 
modes leads to exponential decay of coherence for 
intermediate times. 

• In Section IVII we conclude. In light of the results 
we obtain in this paper, we critically examine the 
degree in which quantum memories can be encoded 
using Majorana fermions when these are in contact 
with a dynamical environment. We show that the 
coherence of the Majorana mode is controlled by 
the coherence of the bath it interacts with. 

• In Appendix |X] we compute exact dressed zero 
modes for static quadratic Hamiltonians, which we 
use to verify the validity of our results in Sec- 
tion [IVl In Appendix [B] we present a rather tech- 
nical calculation of a Majorana mode interacting 
with a fermionic bath with fully quantum mechan- 
ical Gaussian fiuctuations. To leading order we find 
a decay similar to classical computations. In Ap- 
pendix [C] we present various technical calculations, 
used throughout the rest of the text. In particular, 
in Appendix IC II we show that our results are in- 
dependent of coding subspace, in Appendix IC 31 we 
present some technical arguments (which are used 
in Section |V| in favor of weak (negligible) coupling 
of the fiuctuation for the various fermionic modes. 



In the rest of the appendix we derive formulas used 
in the main text. 



III. DYNAMICS 

We begin with a study of the general properties of the 
dynamics of a system of Majorana modes. We will fo- 
cus on a computation of correlators involving Majorana 
operators. This will allow us to study the coherence prop- 
erties of a topological quantum memory which is based 
on qubits made up of localized zero energy modes. In 
this Section we will adhere to very general Hamiltoni- 
ans and we will study only properties that are essentially 
independent of the form of this Hamiltonian. This will 
set us up for studies of specific types of Hamiltonians in 
Section IIVI From the outset, we would like to specify 
the initial conditions or equivalently the density matrix 
when the system is initialized at i = 0. We will assume 
that initially the density matrix factorizes into a product 
of the form: 



Plot = PMaj ® JJ Po 



(8) 



Here ptot is the density matrix for the entire system, 
while pMaj represents and arbitrary non-equilibrium den- 
sity matrix for the Majorana modes. The pcnv; are ar- 
bitrary, not-necessarily equilibrium, density matrices for 
the environments of the individual Majorana modes. No 
specific "ensemble" is assumed. This form is a reasonable, 
consistent assumption for the initial states of system plus 
bath, particularly so, as many experimental methods of 
initialization produce such states. 

For our qubit memory persistence between times ii 
and ^2 is captured by the two-time correlators such 
as (cr^ (ii) cr^ (^2)). We note that, because the initial, 
t = 0, state breaks time-translation invariance, gener- 
ically these correlators are functions of both ii and 
^2. Here we shall focus specifically on correlations, like 
(cr^ (0) (7^ (T)), between the state prepared at i = and 
the state at a later time time i = T which characterize the 
degree to which the information encoded in the qubit at 
the initial time survives interaction with the bath when 
it is retrieved at a later time T. 

The key results of this section, which are used repeat- 
edly later in the text, may be summarized by saying that 
even though the factorization form given in Eq. (|5]) does 
not survive Hamiltonian evolution the expectation values 
of various correlators like (ct^ (0) cr^ (T)) or equivalently 
products of Majorana fermions, to be defined precisely in 
Eqs. ([T^ and P^ below, do factorize into products of 
expectation values for individual Majorana modes. This 
factorization survives for arbitrary times. 



A. General ideas 

We will consider a set of Majorana modes each inter- 
acting with its own fermionic environment, see Eq. ([S]). 
We will see that there is decoherence even without di- 
rect interactions between different Majorana modes or 
between their respective environments. One can show 
that, in the limit when the spatial separation between 
the Majorana modes is large, the case when multiple Ma- 
jorana modes interact with a common fermionic bath re- 
duces to the case of uncorrelated non-interacting baths 
(see Appendix IC2[) . The Hamiltonian pertinent to each 
mode may be written as: 

1=1 

Here Bi^^ are some bosonic modes and a ~ {1, 2, ...} la- 
bels the Majorana modes. The total Hamiltonian is given 
by _ff = J^a^a- We will be interested in correlators 
of the form (7^1 7^2 .. . Ja^ 7qi (ti) 7^2 (^2) • • • 7ak i^k))- 
Here all operators are in the Heisenberg picture, and 
ja {t) is given by 

7„ (i) = (f e' fo H. (^)<:;^) ^^ (j e-» I* H^ (r)d.^ 

(10) 
where T and T stand for time-ordered and anti-time- 
ordered products, respectively. Notice that 7q, (t) = 
7J, (t) at all times. 

Now, by Taylor expanding the time-ordered and anti- 
time-ordered exponentials in Eq. (jlOp . taking various 
commutators, grouping terms and using the fact that 



1, we may write that 



-ra{t) = ic.B^{t) + T^{t), 



(11) 



with Ba (t) and Ta {t) having no factors of 7q,. Because 
"/a (t) must be fermionic (this can be seen from the fact 
that the Hamiltonian and all its powers are bosonic) 
we may deduce that Ba (t) and J^a (t) are, respectively, 
bosonic and fermionic operators. By the conservation of 
fermion parity we know that the expectation value of any 
operator {Ta {t)) ~ 0. Finally, because ja (t) is Hermi- 
tian, it also follows from the properties above that Ba (t) 
and Ta (t) are Hermitian as well. 
Now, it follows that 

(7a7a(i)) = {Ba(t))+{-faJ'c.{t)) 

= {Ba (t)) + (ja) {^Fa (t)) 

= {B. it)) , (12) 

where we used going from the first to the second line 
of Eq. p^ that the environments and the Majorana 
states are initially disentangled so expectation values 
factorize. Note that this comes about because in the 
Heisenberg picture the expectation values for operators 



are taken with respect to the initial state, at i ~ 0. 
For the third Une we have used that the expectation 
value of any fermionic operator (Jq. (t)) should be zero. 



J 



Note that because Ba (i) is Hermitian this implies that 

(7a 7a (t)) e K. 



The following factorization formula can be similarly showed: 



(7a,...7a.7a,(tl)...7a.(i..)) = {-if'-'^/' n(^".(^^)) 

= i-lf'-'^^' n(7«. 7a, (t,)) 



(13) 



for distinct aj, j = 1, . . . ,k. To show this expression, one uses Eq. (|lip and again that the expectation values are 
computed with respect to the initial density matrix given in Eq. ([5]) which has the property that the environments are 
uncorrelated with each other and with the initial Majorana states. We see that this factorization formula is independent 
of the initial state of the density matrix of the bath. As such our formalism captures highly non-equilibrium initial 
conditions. 



B. Qubit memory correlations 



The degree of persistence of memories assembled using Majorana fermions can be quantified by the correlation 
between the qubit state, encoded as in Eq. ([3]), at two times 0, T: 



(a'~ (0) a-- (T)) = -(7i(0)72(0)7i(T)72(T)> 

= (7i(0)7i(T))x(72(0)72(T)) 



(14) 



Notice that the factorization implies that, even though the qubit is defined non- locally using two spatially separated 
Majorana fermions, the decay of the memory is controlled by the product of the two-time correlations of the two 
separate Majorana modes. In particular, the decoherence rate is independent of the initial state of the quantum 
memory (that is correlators of the form (7172) do not enter the result). 



r 



Thus in the case of uncoupled well separated Majorana 
modes each interacting with its own environment the task 
of determining the persistence of topological quantum 
memories based on Majorana fermions is reduced to the 
calculation of the coherences (70, (0) 7q, (T)) in the pres- 
ence of different fermionic environments. We carry out 
this program henceforth. 



use we will study both static and time dependent Hamil- 
tonians. Based on the discussion given in Section Hill for 
the purposes of computing coherences it will be sufRcient 
to focus on a single Majorana mode. As such we will 
drop the subscript a, see Eq. ^, henceforth. 
A. General Observations 



IV. KELDYSH CALCULATION OF 
COHERENCE 

We now proceed to describe the technical details asso- 
ciated with studying dynamics. For generality and later 



J 



We will convert the computation of the Majorana cor- 
relations into a Keldysh calculation carried out using only 
the bosons and regular complex fermions inside the reser- 
voir. (For a review of standard Keldysh techniques see 
e.g. l39l - [4ll .) We will calculate the following correlator: 



(7(0)7(T)) = (7 {fe+'^^omr)dr^^ ^ ^^ ^_.,^^ Hir)dr\^\^ 



(15) 



Here the expectation value is taken relative to the density matrix po at r = while T and T stand for time ordering 
and time antiordering respectively. To make the computations tractable we will assume that po = Pthorm ® PMaj ■ Here 



PMaj is any initial density matrix acting on the subspace of the Majorana modes while pthcrm is the thermal density 
matrix for the regular fermion modes. 

To compute the correlator in Eq. (|15p . we will use Eq. ([9]) and work in the interaction picture with respect to the 

rest of the Hamiltonian iJ^^'^ / {cj , Cj , Bj , _Bj } 1 . We will expand the ordered exponentials in powers of H and collect 

and contract all the 7s to eliminate them. In what follows will show that 



(7 (0) 7 (T)) = ( (f e fo o(r) rfr ^ (^^ ^ /- o(.) d.^ ^ ^ ^^^ 



,-E„/o''0(r'')dT 



(16) 



where O (r) = "^i^i ( Bi (r) c^ (t) — bJ (r) c] (r) J , and % stands for the Keldysh ordering that combines the forward 
and backward propagation, and the index a= t, b labels the two pieces (forward and backward) of the ordered 
product. (Notice though that the operator O (r) in the exponential comes with the same sign in the T and T 
products.) 
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Figure 2: The Keldysh contour determining the coherence of the Majorana zero mode. We consider 2K insertions of our inter- 



action Hamiltonian i^iy'^^^i ( BiCi — BJcl I into the Keldysh contour with ± referring to the forward in time and backwards 
in time branches. Several interaction insertions are shown by dashed lines. To convert this contour to a "regular" Keldysh 
calculation we commute the Majorana modes (7 terms) including the one at r = T till they are all located at r = as shown. 
In the text we describe how to compute commutators appropriately. 



Below we give the essential arguments needed to derive Eq. (J16p . To carry out this program, let us introduce a 

short-hand notation H = j J2i=i ( ^i'^i ^ ^iH ) = 7 (-^'^ ^ B^^c^^) . Now expand Eq. ([T5|) in powers of H, and focus 

on the term with A^b + -^t insertions, with N\^ from the expansion of the 7~-ordered exponential and N^ from that 
of the T-ordered exponential. By fermion parity conservation and using our assumption that the system-bath initial 
density matrix is factorized we know that A^b + Nt — 2K is even. The insertions of our interaction Hamiltonian are 
of the form 



{7} [^J (Be - Btct) it\)] ■■■[^J (Be - B^e^) (^J] {7} [-^7 (Be - B^ct) (t*)] • • • [-Z7 (Be - Btet) (4j] 



bottom insertions 



top niscrtions 



(17) 



We show in curly brackets {7} the modes at r = T and 
at r = 0, to help single them out for constructing the ar- 
gument below. Our strategy to convert this calculation 
to a "regular" Keldysh calculation will be to move the 
Majorana modes (7 terms), including the {7} at t = T, 
by taking appropriate commutators, till they are all at 
the left hand side, adjacent to the {7} inserted at r = 0. 
We will move along the contour ordering direction (see 
Fig. [2]). We will then use the relation 7^^+^ = l to 
eliminate these modes altogether. All that remains is a 
computation of the commutators. Because of the form 
of the Hamiltonian, computing commutators is equiva- 
lent to computing an overall sign for the term in the 
expansion. By noting that the Hamiltonian is bosonic 
we obtain that the overall sign is only due to the anti- 



commutation of the 7's with the q and c] inside the 
(Be — B'l'e'f) terms. We shall move each 7 mode to the 
very left in two steps: we first move the mode at r = T 
to the very left towards r = 0; then we move all the 
remaining modes there as well. 

In the first part of the procedure is to obtain the con- 
tribution of the Majorana fermion inserted at t = T. We 
note that the number of — 1 signs it picks up depends on 
its position along the contour relative to the other modes 
it picks up one —1 sign for very mode it passes so there 
is an overall sign of (—1) ''. 

Now for the rest working from left to right, the first 
Majorana mode that needs to be moved picks up no —1 
signs as it does not pass over a (Be — B'l'e''') term, but 
the second picks up one —1 sign as it passes over one 



such term. Similarly, the third picks up two (—1) signs, 2K ~ 1 factors of —1. The product of these factors yields 
and so forth. Finally the 2Kih Majorana mode (last (^_i^k(2K-i) ^ (_i)^' = (_ 
to be moved, sitting all the way to the right) picks up 



■\N,+N^ 



Thus eliminating the 7's in Eq. pT|) leads to an overall sign (— i)'^t,+^' x (— l)^i», which then allows us to replace 
terms of the form Eq. ([T7|) by 

[- (Be - Btct) (t?)] . . . [- (Be - Btct) (4j] [- (Be - Btct) {t\)\ • • • [- (Be - Btct) (4j] . 
^^ . '"^ . ' (18) 

bottom insertions top insertions 



These are precisely the terms that appear in the series expansion of Eq. (|16p. and therefore we can continue the 
calculation utilizing this expression. We should point out that for complex fermions coming from Majorana insertion 
Tc corresponds to literal ordering on the Keldysh contour, without any fermionic minus signs, because the original 
Hamiltonian was bosonic [this can also be seen step-by-step in going from Eq. (|17p to Eq. p8| ]. This fact leads to 
the modified sign for the fermionic 7^-ordering: 



0{h-t2) 4{h)c,{t2)+e{t2-h) Q(t2)cJ(ti) ,ii,i2 0ntop 



nU+ \n a \ - J 4 (*i) '^^ (^2) , ii on bottom, ia on top , . 

Ci [12} c\ (ii) , ti on top, 12 on bottom 

e{t2-ti) c\{ti)c^{t2) + e{ti-t2) c, (t2)4(ti) , ii,i2 on bottom. 



Now, we turn our attention to the computation of Eq. (J16p . We do so in steps, computing the expectation values 
by first tracing the fermions (ci , c\) and then subsequently tracing the bosonic degrees of freedom. Even in the case 
where there are interactions for the fermions, we can still treat the theory as quadratic in the fermions and include 
the interactions (with photons or phonons) as a coupling of the fermionic bilinears with the mediating bosons, which 
we label by (j>- Alternatively, we may think of the fields fields as Hubbard-Stratonovich decoupling fields^. 

We can thus write 

^ ^^^[\Y. j^^ jy^2 (r4(Bc-Btct)(rn {^c-^^c^){tI)])\ . (20) 

We remind the reader that all functional integrals are along the Keldysh contour. The action 5^ is that of the 
interaction mediator field and contains the dressing from the integration of the fermions, which are integrated out 
first as explained above. The normalization Z is 



J{t{v^,,VBI\ e'^BlB^Bl] ll^Jl^^^^^t 



Z = j MpB.PBl e"^Bii.„i.ii j M|p<^^p0T e*-^*'^"^"!. (21) 

This procedure works because it possible to calculate partition functions. Green's functions, integrate fields out etc. 
along any contour, in particular along the Keldysh contour as used here. We then express the fermionic correlators 
in terms of their Green's function, 

(rc[(Bc-Btct)(rf) (Bc-Btct)(r|)])^^^.,. = -B, (rf) S/ (r|) (To [c. (rf) c/ (r|)] ) 

-BHr'i)Bj{r!^) {%[cHri) c, {r'2)]) 
^ -B(rnG^^e(^f,^2) Bt(r|) 

-Bt(rf) Gt,^{T^,,T'2) B(r|), (22) 

where the Gp ^ {jii'''2) ^'^'^ '-'f ii ('''i'''''2) ^''^j respectively, the electron and hole fermionic Green's function, and we 
have used the fact that the bosonic fields B , B^' can be treated as c-numbers as they are inside the bosonic path 
integral. As stated previously Gp ^ (r", r^) and Gp j^ (r", r^) are slightly unusual Green's functions, with no fermionic 



minus signs (only plus signs), as shown in Eq. ([T^ . Let us define D^ (tJ^jtI) — G 



can then write 



'F,h 



K,^!) 






'2' '!)■: 



(7(0)7(T))=Z-i |fn^B,PBtj e^-SelB^Bt] y/jlp^^p^tj ^ 
^ «^P (-^E /V^^drl Bt(rn i^^lr^rl) B (rl) J 



i-Sf^a^l 



(23) 



We remark that the expression in Eq. (j23p was derived without any approximations. It holds for interacting electrons 
as well, as long as the interactions are included via an external bosonic field denoted by (p above. Furthermore 
we would like to note that though it is not used anywhere in this paper, but a similar path integral formulation 
using Grassmann variables may be done without any decoupling fields, for regular quartic ~ ^^ (x) ^l/^ {x) '^ (x) ^P (x) 
fermionic interactions. A systematic Keldysh diagrammatic perturbation theory may be derived from it. 

For future use we note that to compute the coherence of a Majorana mode it is often enough to compute the four 
diagrams shown in Fig. ([3]). Following Eq. (j23p . their sum may be explicitly written as: 



V[T) ^Y. I d^° / '^^2 B^(rr) Dt{Tt,r, 

a,b -^0 -^0 



BK 




T (b,^ {tI) B, (rl)) X [e (n - T-2) (4 (n) c, (r2)) + 6 {r^ - n) (c, (r^) cj (n))) 

+ f [b,^ (ri>) B, {r^)) X (^ {r, n) (cj (n) c, (r^)) + 6 {n - r,) (c, (r^) ^ (n))) (24) 

+ (B,t (r,t) B, (r^) (cj (n) c, (r,)) + B, {rl) B\ (r^^) (c, (n) cj (r^)))] } 

Here T, T refer to time ordering and time anti-ordering operators. This form places the time ordering or antiordering 

terms (T ( -Bj (rj;) Bi (tj) 1 ) with the appropriate fermion correlators so it can be used directly in calculations without 

having to use a path integral. The factor of two going from the first to the second line comes from a symmetry ri -f-?- T2 
(which also allowed us to simplify Eq. ([M)) above to contain six rather then twelve terms). Because of exponentiation 
of disconnected diagrams, if we can safely ignore higher order correlations among the i?i's, we may write that: 



(7(0)7(r)) = 



-p-i(V(T)) 



(25) 



A quick way to derive the extra factor of ^ in Eq. (|25p above is by noting that it is a symmetry factor associated 
with the ability to permute the two Majorana insertions without changing the diagram [alternatively we can do a 
combinatorial check, or use Eq. ([?5)l ]. 

Let us illustrate with a few simple examples how one can use the expression for the Majorana correlations 
(7 (0) 7 (T)) in Eq. (P5|) to calculate the the decay rates of topological memories. We then deploy this expression 
in detailed studies for fluctuating Hamiltonians in Section FVl 




Figure 3: The four diagrams relevant to calculating V (T) in the main text. We need to sum over four possible orderings of the 
Majorana insertions on the Keldysh contour. The value is given by a sum of terms like Bi^ (r^) _B, (12) {c\ (ri) Ci (T2) ). 



B. Simple examples 



Let us consider simple cases where the Bi are simply constants Ti, switched on at r = 0. In this case the expression 
in Eq. (P5)) simplifies to 



^ „ I Jo Jo 



X exp 



(26) 



where D 



<2) 



p(2) 



r^rl) 



7t(2) 



-(2) 



■(2) 



T2,Ti^, with Gp 1, and Gp g exact 2-point electron and hole Keldysh 
propagators, including the effects of interactions. To be explicit at this level of approximation our formalism handles 
all the dynamics of the (f>a fields but treats fermionic interactions to quadratic order. The . . . stand for terms of order 

©(r^) that involve the 4-point Green's functions G . We shall not do so in this paper, but by including these O(r^) 
and higher terms it is possible to handle all fermionic interactions as well. 

Taking into account all the four cases in the sum over top and bottom insertions ^^ j^, one can write 



a, 6 



T 

drf 



dr* 



j(2) 

■'f 



rtD^K, 



E 



T 

dn 



T 
dT2 



rn({cl(ri),c,(T2)}) r 



r 



(27) 



We now consider a case where this formula will be par- 
ticularly useful. We Consider the case when the bath is 
described by the Hamiltonian 



H 



N 

7E(r»^' 






In this case we have 



N 

E 



^7'C,; Ci 



({cI('ri)'Cj(T2)}) =(5y e 



-iei(Ti-r2) 



(28) 



(29) 



with ti the energy of mode i. It follows by substitution 
in Eq. ^^ and then in Eq. (l23|) that 



or 



(7(0)7(T)) 



(7(0)7(T)) 



.-2EJr,|^|/„^dre- 



.|2 



-4E, 



[l-cos(eir)] 



(30) 



(31) 



If the bath has energy eigenenergies e^ away from zero 
energy (i.e., there is a gap e < \ti\), we may drop the 
oscillating terms in the limit of T S> 1/e, so we can write 



(7(0)7(T))«e 



-4E. 



ITjl^ 



T>1/? 



(32) 



In this case, the Majorana memory decays to T indepen- 
dent plateaus at large times. Thus, as long as the sum 
y^ —^ converges, the memory is retained to a finite ex- 
tent. This result is confirmed by a time-independent re- 
diagonalization in the presence of the F^, which is shown 



explicitly in Appendix [^ where a new exact zero mode is 
calculated. Here we simply note that the finite depletion 
found in this case is a simple consequence of the fact that 
the modes change once the coupling is switched on. Also, 

IP |2 

we compute the sum ^ —^ j and find it to be finite, for 
a specific tunneling model in Appendix IC 4dl 



V. FLUCTUATING HAMILTONIANS 

So far we have studied static Hamiltonians. To gain 
further insight it is interesting to extend our results to 
fluctuating couplings (which may come from time depen- 
dent classical fluctuations or from quantum dynamics). 
We shall focus on three cases, in all three the fermionic 
action is quadratic. In the first case we study we con- 
sider the situation when the Bi are simply replaced by 
classical variables Fj, like we did in Sec. IIVBI but now 
they depend on time. The second case is that when the 
energies ti of the electrons in the bath fiuctuate in time, 
because of environmental fluctuations. The third case is 
a generalization of the first one, where we treat the Bi 
quantum mechanically with their fluctuations governed 
by a quadratic action. We treat the first two cases here, 
and the third, more technical one, in Appendix [Bl 

In the first two cases, one can generalize the expression 
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(a) 



G(Ti,T2) ^ 





Figure 4: The two time correlators of the tunneling amplitude 
Gr (ti, T2) = (r* (ti) r (t2)). a) The shaded region represents 
the actual area of integration for Eq. (|40|) . The darker stripe 
represents the area of large values for the correlator. This rep- 
resents strong correlations in the tunneling amplitudes. From 
this we see that the majority of the integrals appearing in Eq. 
(|40|l come from times when ti = T2. b) A simplified integra- 
tion area. The darkly shaded area of large correlators does 
not change significantly. As such geometrically we see that 
this should not change the values of the various correlation 
functions we are studying. From this it is particularly easy 
to derive the estimates used in Eq. (|43p . in particular the 
linear in T scaling can now be derived by simply changing 
co-ordinates in the integral in Eq. (|40[l. 



in Eq. (|30p simply by taking Ti — > Ti{T) or Ci — > ^ilr): 
(7(0)7(r)) = c-25:.|/o-^rr.Me-/o--.<.|^ 

i 

and then average over statistical fluctuations of the F^ (t) 
and ^iir). 

The computation of the Majorana correlations can be 
greatly simplified as follows. Notice that, for each mode 
I, the argument in the exponential in Eq. (|33|) can be 
viewed as the magnitude square of the position Zi of a 



particle moving in two-dimensions, or alternatively the 
modulus square of a complex number Z^ moving on the 
plane: 

Z,(T) = x/2 / dTT,(T) e-'fo^t^'it) , (34) 

Jo 



with 



{liO)liT))=l[e-\^^\\ 



(35) 



Below we will argue both in the cases of fluctuating am- 
plitudes Ti (r) and energies e^ (t) that the probability dis- 
tribution for the "position" Zi is Gaussian: 



PiZ.) 



1 



27raKT) 



exp 



1 \Z^\ 



2^nT) 



(36) 



with ffi (T) the time-dependent width of the distribution, 
which we will compute below for each case. With this 
Gaussian distribution for the Zi, we can compute the 
average Majorana correlation, 



(7(0)7(T)) = Wj d^Z,P{Z,)e-\^^\ 

i 



exp 



-2E^'(T) 



(37) 



In the last step we assumed that there are many modes 
in the fermionic bath, each making a small contribution 
(or order inverse volume) so we may re-exponentiate the 
product. The examples below are studied using this ex- 
pression. 



A. Fluctuating amplitudes 

The fluctuations of the F^ are assumed to be Gaussian 
distributed according to 



T{{T,{T),Tt{T)])=N-'e 



■5 /r=^ dr^ jr^ dr2 r*(Ti) G-/(ri,T2) r.(T2) 



(38) 



Let us show that the distribution of the P{Zi) is Gaussian, and relate (Ji{T) to the fluctuations of the F^. That the 
distribution T{Zi) should be Gaussian is not surprising since at long times the particle is diffusing. We can write for 
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the characteristic function distribution (Fourier transform of the probabihty distribution P{Zi))\ 

^g-»ifc*x/2/„'rdrr.(r)e-".- ^ ^-ji fex/2 /^ dr T* (r) £+"•- 

Therefore, the distribution P izA is Gaussian, with a variance given by 

af(T) = 2 /dn /dT2 e-'^-^i Gr,(Ti,T2) e+'^^-^^ . 
Jo Jo 



(39) 



(40) 



If the noise correlations are invariant under time-translation, then Gri(Ti,T2) = Gri(Ti — T2). We can expand these 
correlations in frequency domain, Gri(Ti — T2) = /^ dw Gri(w) e~*"''^i~'^^^ 

We proceed to compute erf (T) in Eq. (PD|) for two distinct cases of low and of high frequency noise. 



Case I: Low-frequency noise 

In this case, we shall assume that all frequencies w 
for which Gr^ (w) has significant weight fall below the 
fermionic energies e^. It the follows that 



^?(T) 



, /• ^ 1 - cos e, + uj)T] ^ 






(41) 



We thus arrive at a correlation decay, for the Majorana 
modes, of the form 



(7(0)7(T)) 



exp 



-4^4 / duGvA^) 



(42) 



The coefficient in the exponent depends on the spec- 
tral weight of th e noise. From Parceval's theorem, 
/_ dbj Gri(w) = |ri(i)P, so the prefactor depends on 
the intensity of fluctuations of the couplings Ti[t) in time. 
When the fluctuations are large, for example when the 
Tiit) are tied to thermally induced vibrations in two di- 
mensional systems, there is large decoherence. 

We remark that even in the cases when erf (T — > 00) 
is bounded, the value may be rather large, and the Ma- 
jorana correlation is exponential in this value. Therefore 
keeping the error to within reasonable bounds for quan- 
tum error correction to be applicable can be a tall order. 
In this sense, the Majorana qubit is not necessarily any 
more robust than other proposed qubit platforms. 



Case II: High-frequency noise 

In this case we compute erf (T) assuming that the cor- 
relations Gpi (ti — T2) decay in time, so one can break the 
T1.2 integrals into center of mass: (ri +T2)/2 and relative 
coordinates ti — T2 integrals, and in the limit of large T 
one has 



^f(T) 



T largo 



> 2T Gr(eO , 



(43) 



where Gr(ei) is the Fourier transform of Gr(T) at fre- 
quency Ci. We further clarify this in Fig. Q. 

We thus arrive at a correlation decay, for the Majorana 
modes, of the form 



(7(0)7(T)) « exp 



-4T^Gr(e. 



(44) 



Notice that this expression has meaning only if the 
Gr(a;) has spectral weight above the gap e. If not, one has 
to treat the problem in the low frequency limit discussed 
above. 



1. Non zero expectation values 

One can generalize this result for when the F^ fluctu- 
ations are centered around a non-zero value F^. In this 
case. 



PiZ,) 



1 



2^af(T) 



exp 



l\Z^ 
2 



Zf(T)p 



^f(T) 



(45) 



where 



Z°(T) = ^/2^^ / dr 



V2*r" 



-iciT 



Ci 



which lead to 



(7(0)7(T)) = nf 



|z"(T)r 
+ 2a,(T) 



(46) 



(47) 



Notice that we recover the static result Eq. (|31ll of the 
previous section if there is no disorder [cri(T) = 0]. In- 
deed we see that (7 (0) 7 (T)) = ]J. e-l^°(T)|^_ 

In the particular case of high-frequency noise (non-zero 
Gri(ei))i o'^s obtains in the large T limit one obtains 



(7(0)7(T)) 



T large 



[][l + 4TGr(e,;)' 



(48) 
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which agrees with the case where the fluctuations are 
centered around zero shown in Eq. (|44p . 



2. Cross correlations of fluctuations 

We would now like to extend our model to include cross 
correlations of fluctuations between the modes. Once 
again we focus on a Hamiltonian of the form HMcan = 

7EjIi (rjCi - r*cjj + X^j^i £^c]c^■ Here 7 is a single 

Majorana mode and c^ , cj are regular fermion creation 
and annihilation operators. In our model we will allow for 
Gaussian classical dynamics for the coupling constants Ti 
with possible cross correlations between the couplings. 
More precisely, we will assume that the probability dis- 
tribution of couplings may be written as: 



P ({ri(r), r* (r)}) ^Z-' J Jv {T* (r) , T, (r)} exp I -^ J^^ J^^ dndr, ^ G^j (n, r^) T* (n) T, (r^) j (49) 

Next we introduce the Z = (Zi, Zn) e C^ with Z, (T) = V2 J^ dr T,{t) e-'^" '**'*(*l With this notation we 

may write that: 



-z^z 



(7(0)7(r))=e' 
Which is just a rewriting of Eq. ((35|) . Next following Eq. ((39| we may write that: 

P(/C) = j d^Z^ j d^Z2 j d^Z^.... j d^ZN P{Z)e^i^'^+^'^) 



(50) 



Xe-*3 V^E. 'C.* io ^^ r,(r) e — .- ^ g-ji 721:. K., /.'^ /.'^ dr r*(r) e+" 



exp -i X 2 X V/C*/Cj /dn /dr; 



•2 e""'^iG,j(Ti,T2)e 



+*€ j T2 



(51) 



r 



From this equation we see that the distribution P I Z 
is a Gaussian with a covariance matrix cr (T) given by: 

/.T ,.T 



probability distribution given by Eq. (|49p . In this case, 
we obtain: 



(Tij (T) = 2 [dn fdT2 e-"-"i G„(ri,T2) e 
Jo Jo 

Gombining and simplifying we may write that: 

(7(0)7(T))= ^ 



+ii^j'r2 



(52) 



(53) 



(7(0)7(T)) = 



exp (-f(5 (T) (I + 2a (T))-^ Zq (T) 



dct(I + 2cr(T)) 



dct (I + 2cr (T)) 

rtrix (ly = %). 
the couplings hav 
expectation value, F^ = F? -I- iJF^, with the STi having a Zo^i (T) = yJziT^ - — 7— — ^- 



(54) 



Here, similarly to Section IV A H we have intro- 
generalize to the case where the couplings''have a non-zero duced the vector Zq whose i'th component is given by: 



Here I is the identity matrix (ly = (5.^). We can also 
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B. Fluctuating energies 

Let us consider the case where the energies un- 
dergo Gaussian fluctuations in time, around some av- 
erage value: ei(r) = e; -f- <5ei(r) with {5ei (ti) Sei (T2)) = 

J^dt Se,{t). If 



Gi{Ti, T2). Let ip{t) 



the Sei(T) are short-time correlated the quantity 



[ipin) - (p{T2)] = G^(ri - T2) will grow hnearly in |ti - 
T2I. We note that the phases (pi (r) execute random walks 
in this case. 



The magnitude square of the "position" of the Zi has 
average 



|Z,(T)|2 =2 |r,p fdT+ j dT_ e+"'^+ e+*[¥'(^+)-'^(^-)l e~*''^- 
Jo Jo 



= 2 |r,| 



+ 

Jo 



dT, / dT_ e+"'^+ e"^G'(^+-^-) e"**^'^- 



(55) 



The calculation of higher moments is quite similar if the term 6*^'^'-'^+ '^^' confines the two times to be close to each 
other. 





rjni-p |2n 


/ dr+ . . 


■ r^r^ 


/drr . 


■ /dr,7 




|^,(T)|^" = 


e 






Jo 


Jo 


Jo 


Jo 








/■T 


/■T 


/•T 


fT 






= 2"|r,|2" 


/dr+ .. 


./dr+ 


/ drr . 


■ / dr,7 


e 






Jo 


Jo 


Jo 


Jo 








^j "^j X e *^^ ^j "^j 



X exp 



1 r^ 



dui.. / du„ / dvi.. / dwn ^ {G {ui, Uj) + G (w^, Vj) - G {ui, Vj) - G («», Uj)} 



2"|r,f" n\ 



= n\ (|Z,(T)|^)" 



dr+ / dr_ e+^'-^+ e-5'=^»'(^+-^-' e 



(56) 



For the second equality we have used the fact that the process is Gaussian. In this way we mapped the problem to the 
partition function of a two species Coulomb like gas. Then in the fourth line we have used a dipole approximation for 
the partition function. We note that this is consistent with the confining assumption as / ^ / ^ dudvG (u, w) oc |ri — r2| 
so that we have a confining linear potential between oppositely charged particles of our Coulomb gas. 

We now claim that Zi will execute diffusion because of the random phases. Indeed, these correlation functions 

are the moments of a Gaussian distribution with variance \Zi (T)| . This variance can often be computed in the 
high-frequency case (similarly to Section IV A[) and for large T one can approximate 



\z^mr 



T large 



-> 2T |r,P / dr e+^'-^ e'^^^^^^'^ = TO, , 



(57) 



and the probability distribution is given by P (Z^ (T)) = ,A,„, *^^P ( ^ l To^trT ) ■ Repeating the analysis of 
Section IV A[ we get a power law decay (for each mode i) for the coherence of Majorana qubit, with a coefficient that 
is dependent on the Fourier transform of the exponential of the G'^{t) correlation function: 



(7(0)7(T)) = n 



1 + 4T IF,- 



dr e+"'^ e'i^i'-''^ 



exp 



/oo 
dr e+''^' 



iG^M 



(58) 



For G'i{T) ex |r|, the Fourier transform of e ■^ v^'^' will decay as a power law in frequency. We would like to 
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point out that if the ei (r) have a correlation time to = 
ri~^, the short-time behavior of G'^{t) is smoothened, 
and the kink-singularity of at r = disappears, while the 
long-time behavior \t\ remains the same. Using general 
results on Fourier transforms^ we know that the Fourier 
transform of e^^G'vf'^) will decay faster than any power 
of frequency io when uj ^ fl. This indicates a good level 
of protection for systems with large gaps compared to the 
bandwidth of the noise source. 



between ztA^ with time intervals between events that are 
distributed randomly with characteristic frequency fi" . 
The complex number Zi{T) will again perform a random 
walk at long times, which we will confirm by computing 
the moments of |Zi(T)p. Let us start by computing the 
second moment: 



C. Telegraph noise fluctuations of coupling 
amplitudes 

Here we shall study classical telegraphic noise. Our 
model for telegraphic noise will be a ri(T) that switches 



J 



|Zi(T)|2 =2 dT+ dT^ e+«'^+ r,(T+) r,(r_) 



(59) 



Now, |Zi(T)p = 2Af (— l)^fiip=('^-''^+)j where A^flips (t_ , t+ ) is the number of switches between the two times t±. The 
average 






SO we obtain 



|Zi(T)|^ 



T largo 



-^ 2T Af 



4^: 



{2n,)^ + e^ 



(60) 



(61) 



In the appendix we compute the higher moments and show that the distribution of Zi(T) approaches a Gaussian, 
as intuitively expected from the fact that the telegraph noise causes the fictitious particle position to diffuse at times 
larger compared to the switching time. We obtain, similarly to the previous cases discussed above, that 



(7(0)7(T)) = n 



exp 



1 + 2T A? 



-2T J2 ^ 



4f). 1 


-1 




-1 


'' (2f^.)' + 


A 



(62) 



r 



In the last line we assumed that there are many rele- 
vant fluctuating levels each making a small contribution 
so that we are able to re-exponentiate. From this we see 
that due to the effects of telegraph noise the informa- 
tion stored in the Majorana qubit is lost on a time scale 

~ Ttyp/X^i '~^- Here rtyp '-^ 17"^ is the typical switching 
rate for the regular fermion modes. This is an exponen- 
tial decay of Majorana coherence with the rate given by 
a rational function of the the coupling strengths and fre- 
quencies of the switching. This leads to short lifetimes of 



Majorana modes. We would like to note that the power 
law term comes from the instantaneous switching pro- 
cess. For a finite switching speed and as such a smooth 
(r (r) r [v)) the Fourier transform in Eq. (p^ would de- 
cay faster then any rational function of ti for large e^ (as 
compared to the inverse switching time)ii. 
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VI. CONCLUSIONS 

In this work we have studied the stability of qubits con- 
structed from Majorana zero modes, for example using 
an encoding such as a'^ = 17172- The persistence of mem- 
ory can be measured from two-time correlations such as 
{a^ (0) cr^ (T)), which we have shown is independent of 
the particular state of the qubit. We have shown that 
the if the environments coupling to each Majorana mode 
are uncorrelated, then the qubit overlap function factor- 
izes: (a^(O)aMT)) = (7i(0)7i(T))(72(0)72(T)). We 
then analyzed, in detail, the decay of the Majorana two- 
point function (7(0)7(T)), when the Majoranas couple 
via tunneling to fermions in a bath. We considered only 
baths where the fermions had a gapped single particle 
spectrum (gapless baths would trivially destroy coher- 
ence). We considered both cases where the tunneling 
amplitudes were static, and cases where they were dy- 
namical, fluctuating either classically or quantum me- 
chanically, say mediated by a boson bath. 

Static tunnelings are, expectantly, not consequential 
leading to finite decay. Though this serves as a way 
to check our generic formalism. More precisely if the 
fermions in the bath are non-interacting and if the tun- 
nelings are just switched on but then kept constant there- 
after, then the Majorana qubits only experience a finite 
depletion which we checked by explicitly rediagonalizing 
the non-interacting fermionic Hamiltonian with the new 
couplings. This result can be easily interpreted as a finite 
adjustment in the overlap of the qubit before and after 
the basis changes upon switching the tunnelings. 

However, dynamic fluctuations of the tunneling am- 
plitudes can have very serious consequences. Our anal- 
ysis makes it clear that the dephasing of the Majorana 
correlations is tied hand-in-hand to fluctuations (spec- 
tral functions) of both the fermionic bath and the noise. 
In some instances, for example in the case of athermal 
telegraphic noise, fluctuations can destroy the Majorana 
memories, leading to complete decay of coherence at long 
times. We analyzed several types of noise in the bath, 
both classical and quantum. To understand the rate of 



information loss in experimentally relevant systems it is 
important to study various materials, relevant sources of 
noise and in general realistic spectral functions of the 
bath. The formalism here presented forms the basis for 
such analysis. 
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Appendix A: Non interacting systems (quantum 
depletion) 

To have yet another independent check of the results 
presented in the paper we would like to derive results 
similar to Eq. ([5^ in a different way. More precisely 
we will consider a model consisting of a Majorana mode 
interacting via tunneling with non-interacting complex 
fermionic modes. The Hamiltonian of our system will 
be: 



iiyVc 



N 



r.&-r*c 



N 

E 



Eicjc; 



(Al) 



We will first proceed by exactly re-diagonalizing the 
Hamiltonian. By taking commutators of the form 



[-fflVIcan, 7] J [-fflvlcan, Ci] and 



H, 



Mean 5 C, 



we may rewrite 



this Hamiltonian as a matrix acting on the space spanned 

by "^ ^j Ci, cj > (the factor of ^/2 is a normalization 

constant that insures that the matrix representing the 
Hamiltonian is Hermitian in this basis). With respect to 
this basis we may write that: 



J 



Hmc 



( ^f2Y^ ••• 
\f2Y\ ei 



-\/2ri 



V -%/2rjv 
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V2rjv -\/2rt 










ejv 










-ei 










-£2 











(A2) 



We may now diagonalize this matrix by solving for the 
eigenvalues of the system {A^} with corresponding eigen- 
vectors |K = C/k7 + Eill Unad + J2f^i Uf,.N+i4\. 
By direct substitution into the equation HVk — XkV^ we 
see that: 






c/« 



N^ 



Ak — Ci 



(A3) 



C/. 



Here we have ignored the "top line" of Hmchu in Eq. 
(|A2)) . Substituting Eq. ((A3|) into the "top line" of i^Mcan 
we get that: 



A«C/« = ^V2r,[/,,,-^x/2r*(7,,w+, (A4) 

(A5) 






(A.f-(Qf 



We can now obtain eigenvalue equations: 



A«==0, or 1=4^ 



(A.)' - (e^f 



(A6) 



Now substituting Ao = into Eq. (jA4|) we get that: 

JV N 



\uo\'(^+^f:^^ 



(A7) 



From this we see that the overlap of the new zero mode 
with the original mode stays finite (which would lead to 
non-zero coherence for arbitrarily long times) whenever: 






E 



< oo 



(A8) 



This result is similar to Eq. p2p in the main text. This 
condition is true for any finite system. However the over- 
lap of this mode with the original zero energy mode is 



J 
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depleted by a factor of: 



N 



N 



1 + E 1^0.^1' + El^O'^+'l' 



1 + 4E' 



-1/2 



-1/2 



(A9) 



Below in Appendix IC 4dl we will show that this will re- 
main so for mean field like infinite systems. 

Appendix B: Quantum fluctuations 



We would like to extend the previous results, see Sec- 
tion |Vl to the case where the couplings Ti are allowed to 
have quantum fluctuations. That is we will allow for dif- 
ferent fluctuations for the backwards and forwards time 
paths. Once again we will focus on a single Majorana 
mode which may be well described by a Hamiltonian 



of the form HucaniTi, T*) 



7 



Ef.i ( 



T,Ci 



Si=i ^i^iCi- tlere 7 is a single iVlajorana mode and q, c, 
are regular fermion creation and annihilation operators. 
In our model we will allow for Gaussian quantum dynam- 
ics for the coupling constants F^. We will not be able to 
emulate the diffusion equation derivation given in Sec- 
tion ?? but we will provide a brute force resummation of 
the leading order terms contributing to coherence. The 
key difficulty in modifying the approach of Section?? to 
the case of quantum noise is that because of the various 
theta functions, see e.g. Eqs. ^^ & (|B3p . the fermionic 
part of the correlation function cannot be written in a 
factorisable form Gp {ti, T2) 7^ Gip (ti) x G2F (^2) (or 
a sum of such terms) . As such we cannot simply study 
the diffusion of one or several modes, see e.g. Eq. (|M|) . 
but we have to study the diffusion of an infinite number 
of degrees of freedom (which is more difficult). We now 
proceed with the computation, by using Eq. (|24p we may 
write that: 



(7(0)7(T)) 



AA//p{rt,r}cxp(^-iEa,. /oV/o^rf-2^ rt (g^?) (rf,r,^))" r 

X7exp (ifj^ {i^Mcan (F^ (t) , T (r)) } dr) 7 exp {-iT J^ {Huc.n (Pt (r) , F (r)) } dra 



AA//I?{Ft,r}exp(^-iEa,. /oV/o^rfr| Ft [g^ {T^,r, 



j(2) 



ir?,r^) 



(Bl) 
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Here G^^^ = ®, 



Gil ("^1' ^2) G\2 ('^1' ^2) ] TV" = det G^p' and ZJp ' (rj^, r,^) was defined in Eq. mh. We note that 

G*2l(Tl,r2) G^2(n,r2)y'' ^ F VI' 2; H ^_|/ 

Eq. (|25p does not apply as there are correlations between the F's. As such we must compute a functional determinant 
as shown in Eq. (jBip above. We now use the equation: 



dzi...dzndzl...dz*^exp ( — -z^G -^z 



(27r)" det (G) 



(B2) 



Which is true even for an arbitrary (not necessarily Hermitian) matrix G. We will provide an independent proof of 
this result in Appendix [C) Now noting that the determinant of a block diagonal matrix factorizes and writing out 

the form of Dp (rf , T2) say by using Eq 



we can show that: 



h(o)7(Tr'.n*' •+2fa;!;::S!a!;::S! 



G'2i(t-1,T-2) G^2('^l''r2) 



/ e{u-t2){c\{u)c,{t2)) / u^ \ u \\ 



V 



+ e(t2-tl)(c,(t2)ct(tl)) 

/, N t /, A e(t2-tl)(cl(tl)ci(t2)) 



(B3) 



We have inserted the forms of the various matrices explicitly. What remains is to evaluate the functional determinant 
in Eq. (jB3P above. First by conjugating all matrices above with the matrix ^ J (here I stands for the identity 

matrix on [0, T] x [0, T]) we may write that: 



(7(0)7(T))-^ = 
n,dctJl + 2('^"^ ^i 



= nidet<(l + 2( ^j, ^t 






e(t2-tl){(ci(t2)cj(tl)) 

_-(c+(ti)c,(t2))} 

1 ^ u 

G^ Gf j I 



e(ti-t2){(cj(tl)ci(t2)) 

-(c,(t2)cl(tl))} 

Cr{t2)c\{h)) + {c\{h)c^{t2) 




(B4) 



We would like to note the unusual bosonic minus signs in Gf & Gf in Eq. (|B4p above. The rest of this section is 
an evaluation of the determinant in Eq. (jB4p above. Using the identity det (I + M) = exp (^ ^~Tr (M")) we may 
write that 



(7(0)7(T)) = exp E^Tr E,, 



nc: 



'«2fc-lT*2fc i^^2fc:*2fc + l 



(B5) 



Here i^ = lor2 and (i^, ifc+i) 7^ (1, 1). To proceed further we will now evaluate each of the traces (to leading order 
for large T) . As such we need to evaluate integrals of the form: 



j^dT,j^dT2..j^dT2^ { ( [Gf /^/^ (n - r2) X {9 (r2 - n) /O (n - r^) /!)' 



(B6) 



X Cf'^''^ (t2„-1 - T-2„) X (6* (t2„-1 - r2„) /6I (r2„ - T2„-l) /I) j X 
g-...(r2-r3)-K.|r2-r3| ^ ((^ _ 271,) 6 [r^n-l - T2„) / (2n, - 1) (t2„ - T2„_l) /I) 

.... X "e— .(-2,.-ri)-K.|r2„-ri| ^ ((^ _ 2n,) (n - r2„) / (2n, - 1) ^ (t2„ - n) /I)] ) } 

Here for future convenience we have written out the various theta functions involved and for simplicity assumed re- 
laxation time approximation for the fermion Greens functions. The terms AjRjK refer to advanced/retarded/Keldysh 
Green's functions while the various options for the theta functions shown in the brackets correspond to the respective 
green's functions {A/R/K). We now need to evaluate these integrals. As a first step we take advantage of the short 
range of our correlation functions (see Fig. ^) to change range of integration limits for the variables ri, T3, ...T2„_i 
from (0, T) to (—00, 00). We also shift the variables of integration calling Ui = r2i-i — T2i, vi = T2i. Combing all these 
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changes we get that the any term in expansion in Eq. (jB5[) e.g. Eq. (|B6[) may be written as: 

OO />00 /"OO 

dui I du2.... dun X (Gf/-^/^(ui) xe-''^-"i-'='l"il X (6'(-Mi)/6l(ui)/l)j X 



OO ^ — OO 



(B7) 



X (cf/^/^K) xe-"'""-«'l""l X (0(-«„)/^K)/l)) 



X / dui / dv2.... dVn X (6'(W2 - Wl +M2)/6'('yi - W2 - W2) /I) (6'(W1 - «„ +Ml)/6'(w„ - Wl - Ml) /I) 

Jo Jo Jo 

We may further simphfy this expression by noting that all the correlation functions G,- ' ' are dominated by small 
values of u so that we may approximate 9 {v2 — vi + U2) = {v2 — ui) and similarly for other 9 functions. Substituting 
we get that the integrals simplify: 



. J=i 



n/ Gf/^/^'c 



^/"'/^ („^^-^^^^J-^^WA 



^{-u,)/9{u.0/l)}x- 



/ dvi.. I dVn Yl {Q {V] + l - Vj) 19 {Vj - Vj+i) /I) \ 

Jo Jo j^i I 



(B8) 

In Appendix [C] we will further simplify the expression in Eq. (|B8|) above. Here we will merely compute the 
leading order term for the semi classical case where Gf ^ Gf^Gf. In this case a single term (containing only Gf 
contributions) dominates at each order of integration and we may write that: 






(B9) 



Here Gf- (e^ — iKi) is the "Fourier transform'' of the Keldysh Green's function evaluated at energy ti and decay term 
Ki. Combining these results we recover the semiclassical result that: 



(7(0)7(T)) = W 

i 



1 + 2TGf (e, - IK,) 
1 



(BIO) 



1 + 2T f G, (e, - iK,) + G, {-e, + in,) 
exp I -2T^ (g; (e,; - in,) + G, {-e, + in,)^ J 



r 



In the second step we have used a relation between 
Keldysh and time ordered correlation functions and in 
the last step we have assumed that there are many rel- 
evant fermionic modes in the bath so that we can safely 
exponentiate each term. Further corrections to this re- 
sult are given in Appendix [Cl 



Appendix C: Various Tedious Calculations and 
Proofs 



1. Parity eigenvalues (coding subspace) 

In the main text (see Section |l| we presented 
a specific encoding of the Majorana qubit that 
used the even Majorana fermion parity subspace 
for its coding space. Throughout the main text 
we computed expectation values of the form 



(71 (0) 72 (0) 71 (T) 72 (T)) = - (a^ (0) a^ (T)) . We 

claimed that this is a good representation of the fidelity 
of our quantum memory. There could be further 
concern that we are over or under estimating the fidelity 
by including in the expectation value (7(0)7. ..7(T)) 
processes that included final states that do not have an 
even fermion parity^. Here we show that for two time 
correlation functions such processes never contribute 
to this expectation value so no further measurements 
or corrections are needed to adjust for such processes. 
Even though we do not focus on this case in the main 
text we will show that the above statement is not correct 
for multitime correlators. We will also show what 
modifications must be made in the multitime case. 
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a. Two time correlators 



We start by showing that no modifications are nec- 
essary in the two time correlators case (again focusing 
on the four Majorana fermion qubit). Indeed consider 
Y\, and Y\_ projectors into even and odd Majorana 

fermion parity subspaces {X\+ + H- =^ 1' Y\± = Y\± ^"^^ 



Y\,Y\- — 0)- Since the initial state of the Majorana 
qubit has even fermion parity, we may write that: 



(n+'^Mo)n+^MT)n+) 
(ano)n+^MT)n+) 



) 



(Cl) 



In the third step we have used the fact that 
K (0) ' n±] = to get rid of the term 0+ cr^ (0) 11- = 
0. From this we see that we may as well 

project out the odd fermion parity subspace, e.g. 
cr^ (T) — > W^ a^ (T) Y\^ and not worry about errors in- 
volving non-coding subspaces (these errors do not con- 
tribute to expectation values). The same sort of argu- 
ment may be made for any two time correlator of the 
fermion modes and any encoding subspace. Indeed based 
on the form of the previous proof to ensure that the non- 
coding subspace does not contribute to the expectation 
values all we need is a coding system such that the logic 
operators do not take us out of the encoding space (which 
is always the case). So no further corrections are needed 
in this case. 



b. Multi-time correlators 

In the multi time case in order to only consider terms 
within the even fermion parity subspace it is necessary to 
project out the odd fermion parity states explicitly; that 
is convert Oi (T) -> n+ ^i (T) n+- There are still many 
simplifications in the case of three time correlations. In 



this case similarly to what we did above one can check 
that it is only necessary to project out once just before 
the last operator. For example: 

"2^ (717271 (n) 72 (n) (1 + 71727374) 71 {'^2) 72 (T2)) , 

which we can calculate using the methods derived in this 
paper. 



2. Cross Correlations between Majorana baths 

In the bulk of the text we have discussed the case when 
the different baths surrounding the Majorana fermions 
are uncorrelated, or equivalently that interactions be- 
tween modes that couple to different Majorana fermions 
are negligible. In this section we shall discuss the ef- 
fects of such interactions, and indeed argue that they 
may well be neglected in the case of well separated Ma- 
jorana modes: modes whose separation is much greater 
then the scattering length in the bath medium. 

First we begin by arguing that the initial condi- 
tions which we have selected in this paper, of uncor- 
related distant baths, are likely to be highly favorable 
for the coherence of a qubit composed of Majorana 
fermions. Indeed, focusing on two Majorana modes, we 
note that the coherence of the qubit may be expressed 
as (7172 e*^"^ 7i72 e~'-^"^). We now consider two Majo- 
rana modes each interacting with the same fermionic en- 
vironment: in particular we will focus on a shared modes 
/j with energy e, coupling to both 71 and 72 through a 
Hamiltonian of the form i? = 71 Ee (^1 fe - Tf * /j) + 
72 Se (ri /c ~ r|* /I). Here Ff 2 are just complex tun- 
neling amplitudes, for simplicity. Taylor expanding the 
exponentials in the equation above, we obtain non-zero 
contributions to the coherence (the expectation value 
given above) that contain cross terms involving both of 
Tl and r|: 



-2 (7172) pdh l\ J2 ([(^1 Mt^) - rr nih)) , {ri Mt,) - Ff /l(t2))]) 

Jo Jo £ 

= 2(7172) Ail A2 J2 rrr^ ((/.(i2)/i(ii))-(/Kii)/.(i2)))+h.c. 

Jo Jo 



(C3) 



r 



These are the interference terms that do not appear for 
Majorana fermions interacting with separate baths, but 
appear due to a common bath. For short times any non- 
zero terms like those lead to decoherence. Indeed, since it 
is impossible to have higher then unity coherence, these 
terms must contribute negatively to the performance of 
a qubit composed of Majorana fermions. 



However we would like to now argue that this effect 
can easily be avoided in realistic experimental situations 
by simply keeping the Majorana fermions far apart. First 
note that individual / modes that are localized cannot 
have large tunneling overlaps with two distant Majo- 
ranas, so FiFj = 0. Therefore only extended modes 
can contribute to the interference terms. Now, each such 
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mode contains a normalization factor proportional to in- 
verse square root of volume, so individually they con- 
tribute zero in the thermodynamic limit. As such, in 
order to get a non-zero value for the term shown in Eq. 
(jC3p we need to integrate over the contributions of all the 
extended states. To do so first recall Eq. (j4]) or Eq. ()C6P 
below which state that ^l 2 ^ J dr ui^2 {r) x w^ (r). Here 
Ml. 2 is the wavefunction of the Majorana mode while Ug is 
the wavefunction of the mode fe- Assuming a pointlike 
Ml. 2 or dividing the integral into portions of negligible 
extent we may write that T\ 2 cc v^: (ri_2), where ri 2 are 
the locations of the two Majorana modes. In this case, 
we can relate terms entering Eq. (jC3[) to single-particle 
Green's functions for the bath electrons: 

^rrr^(/e(t2)/l(ti)) 

e 

(X^M:(ri)M,(r2) (/e(i2)/l(il)) 
e 

= G(ri,ti;r2,t2) • (C4) 

In a realistic material there are always sources of decor- 
relation, in particular lattice disorder and phonons. It is 
not too difficult to show thali^"— these sources lead to 
an exponential decay of G (ri , ii ; r2 , ^2) in space with 
a characteristic length given by the mean free path of 
the material. The mean free path is directly related to 
phonon and impurity scattering strengths^ir— . Since 
this reasoning indicates an exponential suppression of 
these interference effects with distance, and since it is 
not possible to use these interference effects to enhance 
coherence anyway, we have ignored the possibility of the 
Majorana modes sharing a common bath in the text. 



ent modes of our environment, we mostly focused on the 
case of diagonal fluctuations. Furthermore our results on 
cross-correlations are technical and in practice difficult to 
apply. Here we shall present a partial justification indi- 
cating that diagonal fluctuations are dominant over cross 
correlations. Weak correlations do exist so no "theorem" 
indicating a lack of cross-correlations can be presented. 
We will however present arguments supporting indepen- 
dent correlations in three key cases: when there is a high 
degree of symmetry for the problem, when there is "dis- 
order averaging" of the continuum states and tunnel cou- 
plings have short correlation length, or to leading order 
in perturbation when the fluctuations are weak. 
a. High degree of symmetry 

Many Hamiltonians have a high degree of symmetry. 
For example for a p-wave superconductor with a sin- 
gle vortex supporting a single Majorana mode the vor- 
tex core states have rotational symmetry. Most exter- 
nal Hamiltonians causing fluctuations in the vortex core 
are invariant under this rotational symmetry and as such 
they may be written in block diagonal form with each 
block corresponding to a different eigenstate of the rota- 
tion operator. As such fluctuations corresponding to dif- 
ferent angular momentum eigenstates are decoupled from 
each other (uncorrelated), justifying this assumption in 
this case. More generally fermionic modes corresponding 
to different irreducible representations (diagonal blocks) 
of some fluctuation Hamiltonian have uncorrelated fluc- 
tuations. This in part justifles the assumptions used in 
Section M 



Partial justification of independently fluctuating 
modes. 



In Section |V] we presented some results for the coher- 
ence of a single Majorana mode in the presence of a 
fluctuating environment. While we covered both diag- 
onal fluctuations and cross correlations between differ- 



J 



b. Short correlation length & disorder averaging 

We shall now focus on a particularly simple, but re- 
alistic, model of tunnel couplings between the Majorana 
mode and the regular fermion modes in the supercon- 
ductor. We shall assume point like tunneling with an 
effective coupling that may be written as: 



+ 4 [ (fr {S (r, r) uq (r) v* (r) - S* (r, r) vq (r) u* (r)} 



(C5) 



Here Ui (r) and Vi (r) are the creation and annihilation components of the modes Ci while mq (r) and mq (r) are the 
creation and annihilation components of the Majorana mode and S is a tunneling amplitude. For a similar coupling 
form see e.g. Eqs. (|C25p . & (|4]). From this we see that within our model the coupling functions in Eq. (|28)) is given 
by: 



Fj (r) = / cPr {.=, (r, r) uq ir) Ui (r) - S* (r, r) vq (r) Vi (r)} . 



(C6) 



The correlation function is given by: 
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(r:(ri)r,(r2)) 



(fri (fr2{'B{ri,Ti)^*{r2,T2)uo{ri)voir2)u,{ri)u*{r2) + 
+ S (ri, Ti) S* (r2, T2) Mo (''1) vo (''2) Wj (''i) Wj* (''2)) 
d^^Fin, T2)(|uo(r)|'M,(r)u*(r)) + ^* (n, T2) (|mo W|' w« (r) w* (r))} 

rfV{i^(ri,T2)(|Mo(r)|'W,(r)%) + ^* (n, r2) (|mo (r)|' V, (r) <5,,)} . 



(C7) 



r 



Here we able to simplify our expressions by assuming that 
(S (ri, Ti) S* (r2, Ta))^ ^ F (n, T2) (5 (fi - r^) for some 
F (ti, T2) and that (^ (fi, n) S* (r2, T2)) ^ 0. We have 
also performed a disorder average over the bath states 



Ui (r) Uj (r) 
uum states. 



This averaging works well for contin- 



c. Weak Fluctuations 

In many situations there are many fermionic modes 
responsible for the decoherence of the Majorana mode 
and the coupling to any one mode is quite weak. In 
this case even if the fluctuations between the different 
fermion modes are strongly cross correlated the diagonal 
correlations dominate decoherence. Indeed, to show this 
we first recall the formula for the coherence of a Majo- 
rana correlator given in Section IV A 21 (7 (0) 7 (T)) = 
det~ (I + 2cr (T)). We now simplify this formula. First, 
letting the eigenvalues of cr be {Ai}, we obtain that: 



i 

= exp(-2^\^) =cxp(-2Tr(cr)) 

In the second step we have assumed that many eigen- 
values contribute to the product so we can exponen- 
tiate. From this we see explicitly that in many cases 

I 



with weak fluctuations only diagonal terms of the ma- 
trix a matter. These are one particle terms era (T) = 

2/0^ dri J^ dT2 e-''^'^i Gi(Ti,T2) e+'^'^^ and as such 
are much easier to handle. 



Proofs and clarifications of Eqs. (|B2|). (|B8|) . Si 
([62]) 



a. Eq. (MM 



Here we wish to prove Eq. (jB2p for arbitrary (not nec- 
essarily Hermitian) matrices. As a flrst step we wish to 
prove an analogous expression for real Gaussian integrals. 
More precisely we wish to show that for an arbitrary pos- 
sibly complex n X n matrix M and an integral over R" 
we may write that: 



dxi...dxn cxp 



1 



-x' Mx = 



(2^)' 



1/2 



det (Mdf^ 



1/2 



(C9) 



To prove this we first note that ^ • xiMijXj — 
^'Y^Xi{Mij ~\- Mji)xj. As such we may safely trans- 
form M ^ H-^ + -^^) ■ ^^^^ ^^ ™^y ^^® Takagi's 
decomposition for symmetric matrices^ to write that 
\{M + M^) = UDU^. Where C/ is a unitary matrix 
and 13 is a diagonal one. From this we see that 



dxi...dxn cxp 



1 



(2^) 



^/2 



(27r) 



n/2 



-X Mx - 

2 I {dct{D)f'\\ei{U) (det(i(Af + MT)))'/' 



(CIO) 



The extra factor of dct ([/) comes from the Jacobian of the change of variables. To proceed to the complex case we 
begin by writing z — x + iy, z* ^ x — iy. Then we may write that: 



z^G-^z^ 



-T -T 

X y 



G-i iG-i 
-iG-^ G-^ 



(Cll) 



As such we may write that: 



dzi.. dzndzl.. dz^ cxp [ ^—z*G z 



dx\..dxridy\.. dyn cxp 



-^{^ y 



-zG-i 



zG-i 
G-i 



(27r) 



dot 



G-i 
-^G-l 



iG-i 
G-i 



G-i 

-zG-i 



^G-l 
G-i 



Next we note that: 



1 { G-^ +G 



-IT 



2 V -^G-i - G-i^) G-1 + G-" 



1 i \ / G-i \ / 1 -i 

1 j (^ ^ (G-i - G-i^) G-i^ j I 1 



Since 



det ( Q 1 j = dot f Q M = 1, dct f _, .^_i _ G^i^") G^^"^ ) ^ '^'^^ '^'"^ ^-^ '^^^ ^^ 



-1T\ 



We get that 



dzi...dzndzl...dz*„cxp (^\^G-^^j = (27r)" (dot (G^^) dct (G^^^)) ^''^ = (27r)"dct (G) 
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(C12) 



(C13) 



(CM) 



(C15) 



This reproduces Eq. (jB2 



b. Eq. 



Here we would hke to further simphfy the sums in 
Eqs. (jBSp and (|B5[) as well as obtain more accurate 
estimates. We begin with Eq. (jBSP above. By con- 
sidering the form of the indices in the trace we see 
that we may represent any term in the expansion for 

GfW Gf\y^ 

ken lines with periodic boundary conditions with each 
I 



Tr 



as a set of bro- 



line representing an appropriate Green's function (see 
Fig. ^). In the quasi classical limit the biggest 

contribution comes from the term Tr < (GfGf) > ~ 

T" (^Gf {ci — iKi)) . The last equality may be ob- 
tained by noting that the various terms in Eq. (jB8[) 
factorize. By noting that most of Eq. (jBSP factor- 
izes we may compute the subleading term including 
combinatorial factors in the semiclassical expansion, it 

is f T" (Gf (e - JA.,)) ""' {Cf (e - ^n^) + Gf {e + ik,)^ 
(for n > 1). This term would correspond to diagrams 
(c)-(f) in Fig. (O. As such we obtain that: 



(7(0)7(T)> = Y[exp}Y^^—^T-(Gfie,-^n,) 



^Tl = 



E 



^-^T" (Gf (e. - iK,)) (G^ (e, - ik,) + Gf (e, + in,) 



n 



r 1 + 2TGf (e, - iK,) 
[|cxp(-2TGf (e, -m,) 



exp ( -^T ( Gf (e, - zki) + Gf (e, + in,) 



n 



cxp 



r 



l + 2TGf (e, -iK,). 
Gfie,~iK,) + Gfie,+iKi) ^ 
4Gf {e, - IK,) 



(C16) 



In the final step we have taken the large T limit. As leading order quantum correction, 
such we recover the semiclassical approximation and the 



(a) (1,2) (2,1) (1,2) (2,1) 

(b) (2,1) (1,2) (2,1) (1,2) 

(c) (1,2) (2,2) (2,2) (2,1) 

(d) (2,1) (1,2) (2,2) (2,2) 

(e) (2,2) (2,1) (1,2) (2,2) 

(f) (2,2) (2,2) (2,1) (1,2) 

(g) (2,2) (2,2) (2,2) (2,2) 




Figure 5: In this figure we consider the second order term 
in Eq. ()B8|I above. We picture the seven terms contribut- 



ing to Tr 



Gf \ / G« 



Gt Gf ) \ G^ GK 



with hues con- 



necting indices in the Keldysh matrix, e.g. (1, 2) stands for 
(^(1,2) _ qR gg^gjj entry corresponds to a Green's function. 
The biggest term contains four Keldysh Green's functions 
(pictured last (g)). The six subleading terms are also shown. 



c. Eq. ISM 



We would like to derive Eq. ((62| . As a first step we will 
calculate the n-point correlation function for telegraphic 
noise. We will find that it is short ranged and this will al- 



J 
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low us to calculate the distribution of the "displacement" 
field Zi (T) (see Eq. ([M]) ) within the dipole approxima- 
tion. We will find that the distribution is Gaussian at 
which point Eq. (jHH) will follow. First we motivate the 
dipole approximation used in Section IV CI To do so we 
compute the n-point correlation function for tunneling 
amplitudes acted on by telegraph noise and observe that 
it is exponentially short ranged. That is we extend Eqs. 
([59)) & (pO)) from the main text by showing that for the 
i'th mode, ti < <2 < •■• < tN, and for N even^: 



N 



N 



nr(i,) =Afexp --^(t2,-f2,-i) 



\i=i 



j=i 



(C17) 

To do so we first we recall the result that for tele- 
graph noise the probability of having exactly K flips in 
some set of interval whose total length in L is given by 

■^ (^) exp (-^)—- Now we know that n.^^r (n) = 

±A^ depending on whether an odd or an even number 
of the r (t;) = —A. At this point it is a straightforward 
combinatorial argument to show that: 



{#r in) = -A} = {Ef=i # Flips in [t2j-i,t2j]} 

(mod 2) 

(C18) 
Combing these results we get that: 



tN 



exp (-2i) 

(C19) 
Here L = J27=i (^2j — i2j-i)- As such we obtain the 
result in Eq. (jC17p . Now we wish to calculate 2n point 
function of the displacement field, see Eq. ([M|) . It is 
given by: 



-.N 



Z,(T)|'") = 2^^Jv{T,iT,)}P{T,{T)}J dT,...J dT2„J]exp(^fcie,Tfc)mr(T,)\ 

{2n ( I I 

^^-l)^"-'expK]{(^P,„(,.)*e, + <5)+2(-l)^"f7,} 
1=0 I \j=l 



T Vx 



(C20) 



vJ=i 



2n 

n 



EL=, {^p.,.{k)ie^ + 5 + 2 {-if n^ j \,ti+i Ei=i+i {^p.,dk)^^r + S + 2 i-i)" n,) 



Here {Ti (t)} refers to the space of all path alternating We will derive the second part of this equation separately 

between -l-A^ and — A^ and P {Fj (t)} is the probability of below. The limit: lim5_>.o comes from the fact that some 

( 1, k < n of the denominators may turn to zero without an extra 

such a path, and we have introduced Ui. = < , , . r , f r ai u n j- ^ ^u j. j-u 

^' 1— l,K>n factor ot o. Also we would like to note that there is a 



sum over the permutation group acting on 2n elements: 
P2n which is there to count all the possible ordering of the 
times {ri, ...r2„}. Now consider the formula in Eq. (|C20p 
as a function of (5 € C It is a meromorphic function, and 
it is not too hard to see that it has poles of order at most n 
(this comes directly from the structure of the denomina- 
tors) . On the other hand we know that for i5 close to zero 
i2«\ ^ 22"A2"T2". This is not obvi- 



order T" from Eq. (|C20)) are those 



^ 



P2„(2fc) 



-^ 



S" ' 
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or ones where 



P2„(2fc-i; 



for fc = 1, 2, ...n. From the fact 



thevalueof nZ, (T)l' 
ous from Eq. (|C20p but is obvious from the definition of 
\Zi (T) 1^". As such all the poles in Eq. ((C20)) h ave to can- 
cel. Now, schematically a typical term in Eq. (jC20p may 
be written as a^j^r- (with A S U N). As all the poles 

in S must cancel we may safely replace a^j^ > ct ^, . 

From this we see that for large T to leading order in T; 



Z^{T)\ 



2n 



T". The only terms which contribute to 



that the correlation function g-2ni|ri-T2| jg gj-^ort ranged 
and from the fact that the phase factors in Eq. (jC20[) 
have to cancel pairwise we see that it is good enough to 

evaluate ( | Zi (T) | " \ in the dipole approximation. From 

this we see that /\Z, (T)|^"\ = n! /|Z, (T)|^\". These 

are the moment functions of a complex Gaussian. Re- 
peating the analysis of Section IV Al we get a power law 
decay (for each mode i) for the coherence of Majorana 
qubit, and Eq. (p^ follows. 

Eq. (C2^) : We now wish to derive Eq. ((C20| . By 
considering the form of Eq. (|C17P and the fact that Eq. 
(|C20p has a sum over all permutations of 2n elements we 
see that its enough to derive that: 






X I n^i — I X I n ^ ^ 



(C21) 



To make this formula easier to understand we write it out explicitly in the case when K — 4. 



dTie"i^i / dT2e"^^= / drge"^^^' / dT^e"^''' 



a\ (a\ -|- a-}) {pi\ -I- a2 + 0^3) (ai -f 0:2 + 03 + ^4) 

pQlT 



0L\a2 (02 + 03) («2 + "3 + ^4) 

g(ai+a2)T 

(a.\ -I- 02) 020(3 (03 -I- 04) 
g(Qi+Q2-(-a3)T 

(Ol -1- 02 + O3) (02 + O3) 0304 

g(Q!l+Q2-|-a3-|-Q4)T 

(oi 4- 02 + 03 4- 04) (02 + 03 -I- 04) (03 -I- 04) 04 



(C22) 



We shall derive Eq. (|C2ip by induction: 

T 



Jo 



•^0 /=i [ \j=2E/c=j"fc/ \j=i+iT,k=i+i 

yliK-iLi:U-^^-i) ( n^i — I X I l\ ■ ^ 



All that remains now is to show that: 



i^n— i — +f2 (-1^-' (ri^^ — I X ( n -T^ 

i=i Ej=i "j 1=1 y yj=izZk=] '^k j yj=i+i Efc=(+i 



Ofc 



(C23) 



(C24) 



To see this equality consider the left hand side of Eq. (|C24|) as a function of oi G C This expression is a mero- 
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morphic function C — > C which goes to zero at infinity. 
By inspection, as a function of ai, it has at most simple 
poles. It is straightforward to compute the residues at 
any of these poles and see that they are all zero, that is 
the expression is actually analytic. We can now apply 
Lioville's theorem^ to conclude that the function on the 
left hand side of Eq. (|C24|) is identically zero. 



for say the surface sates formed when an STI is placed in 
proximity to an s-wave superconductor. Furthermore by 
taking the limit of a zero gap s-wave superconductor or 
by ignoring coherence factors we may model insulators 
or metals in contact with p-wave superconductors. We 
shall assume a constant point tunneling contact so that 
the relevant tunneling Hamiltonian may be written as: 



Summation of Eq. ^Al\) for quadratic Hamiltomans 



dPrT (vl/t^ (r) vl/,^^ (r) + vj/t^^ [r) vj/p^ (r)) (C25) 



We will give an approximate calculation of the sum 
(|A7|) for tunneling into a 2-D superconductor. To con- 
sider a simple example we will focus on the case where 
a p-wave superconductor is in close proximity to a 2-D 
s-wave superconductor with the chemical potential of the 
p-wave superconductor set inside the gap of the s-wave 
superconductor. This is a reasonable simplified model 



This form comes from the fact that for a p-wave super- 
conductor the vortex is in one spin species only, say spin 
up. 

We begin with a review of the relevant wavefunctions 
for zero modes of a p-wave superconductor. The eigen- 
values of our Hamiltonian correspond to solutions of the 
following BdG equation: 



J 






\{^{^ 



i{A*(r),p, + ipJ 




(C26) 



Here A (r) ~ exp {iO) A {\r\), with A (|r|) = -y-Aoo for \r\ < S, and A (|r|) = Aoo for |r| > ^ (we have neglected an 
irrelevant overall phase factor). Here ^ is the penetration depth and Aco is the magnitude of the order parameter far 
from the vortex. From previous studie a^^'^" , for rotationally symmetric type II superconducting vortices, we know 
that there is a zero mode for the Hamiltonian given in Eq. (|C26p . It is given by 7 = / dPr {uq (r) ^ ( r ) + VQ{r)'^^ (r)) 
with: 



uo{r) 
Voir) 



-^Jo (kpr) exp {-x (r)) { ^ _\ 



(C27) 



r 



Here kp = ^/2rnJI is the Fermi wavevector, Jq [kpr) is the 
rth Bessel function and x{r) = f^ J^ A (r). Where A (r) 
is the position dependent order parameter. Furthermore 
a good approximate value for the normalization constant 

is given by iV = 0.06 (if) (seeiS). 

Next we will recall the form of the wavefunctions for an 
s-wave superconductor. For s-wave superconductors we 
may write Bogolubov de Gennes equations in the form: 



2m 

A* 




E 



fir) 

gir) 



(C28) 



equation are of the form: 



(+.-) 
,(+.-) 



/ 



(r) 
(r) 



Ai+^-)e'^^Ji{qr) 
Bi+^-)e'"'Jl (qr) 



(C29) 



Here C is a size dependent normalization constant with 
i = ^ (where R is the system radius). Eigenenergies 
and eigenfunctions are now given by: 



{A+, B+) = 
[iA-,B-)^ 



r_ 

2m 



fJ- 



A2 



(cos [6/2) exp (iip) , sin (0/2)) 
(- sin (9/2) exp (iip) , cos {6/2)) 



(C30) 



Here the top component represents creation operators 
for spin up while the bottom component represents an- 
nihilation operators for spin down fermions; Jl and A 
are the chemical potential and the gap of the s-wave su- 
perconductor. Furthermore a similar equation may be 
written with the spins interchanged and A ^- — A. We 
will place the origin of co-ordinates at the center of the 
vortex in the p-wave superconductor. Solutions for this 



Here tan (6) = 



A 



A 
A* 



exp {i2ip) and Ji are the I'th 



Bessel functions. There are completely analogous equa- 
tions for the opposite spin, with appropriate sign and 
phase changes. Using Eq. (jC6p as well as the symmetry 
between the upper and lower component of the solution 
for the zero mode, see Eq. (jC27P and various symmetries 
between the spin species we see that various trig func- 
tions (such as the sine, cosine and exponential appearing 
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in the solution of Eq. (|C29|) above) cancel out. By taking 
the thermodynamic limit we can convert the sum (|A7[) 



into an integral of the form: X]i=i ~^ 



Bvr / dq 
Jo 



V 



(m-a^) + \/(I^-m 



(a^-m) 



Is-M 



N-" 



T 



/ 



drruQ (r) Jq (gr) 



(C31) 

We note that because of rotational invariance only J;=o terms contribute to the sum. Here uq is the upper 
component of the Majorana mode wavefunction (Eq. (JC27P '). We wish to evaluate the integral given in Eq. (|C3ip 
above. We will begin by evaluating J„ drrug (r) Jq (qr). As a first step we will use the approximate relation that: 

Mo (r) = A^exp ( —-^jr^ ) Ja {kpf) (see Eq. (|C27p and discussion that immediately follows). Next we write that: 



drruo (r) Jq [qr) 



N 
2^ 

N 

N 



dx / dyexp ( -T—^r^ ] Jo {kpr) Jo {qr) 



dx / dyexp — 



A 



. d0, I d02CXp(-M(g(0^)+^-(02)'' 

[2tt) Jo 'n V 4^ 



dx 



: exp 



N 



2^2 7_^ yr^^ 

N , fq^\ 



^^ 'kj, + q^ + 2qkFx) 



4A 

4A ■ 



(4 + 9' 



^ /A f kpC , , ,2 

t;— X 1/ — J X exp — —-:- [q — kp) 
27r \ TTqS, * ^ ^ ^ ^ 



4A 



^g,^g-i/cF(8l)-r- / (i5)2e-*9(e2)-r ((^32) 



r 



Here kp(9i) is a vector with magnitude /c^? and direc- 
tion Oi along the x-axis and similarly for q{02)- In the 
second line we have used a representation of the bessel 
function: Jq {qr) = -^ L^ d9e~^^^^^'^ and r is along 
the y-axis. Here /q is a modified Bessel function of ze- 
roth order and in the last step we have used an asymp- 
totic form of the modified Bessel function /q 
. 2^ 



■ exp 



(*)■ 



This asymptotic form fails near 



q ~ where it should be replaced by /o ( ^ ) = 1 + 

3 (^) + ■■■ ^^ ^^ straight forward to check that this 
correction does not effect the final answer see Eq. (jC33[) 
below. Indeed because of the exponential decay we may 
safely approximate: 



drruo (r) Jo {qr) 




M<fcfi 



(C33) 
From this we see that the integral given in Eq. (jC31|) 



above has effectively a finite range of definition and no 
singularities. As such it is clearly finite. Very simi- 
lar arguments may be used to show that the sum (|A7[) 
is bounded for tunneling contact with any gaped mate- 
rial such as an insulator with the chemical potential of 
the p-wave superconductor lying within the gap. Indeed 
quite generically for an itinerant system we may write 
the Hamiltonian a.s H = — -^^-^ + ... which means that the 
eigenvectors of H wee similar to those of an s-wave super- 
conductor so the integrand in Eq. (|C31[) above also has 
exponential decay for large momentum as the solutions of 
H \'i) = i? 14") would behave almost like Bessel functions. 
Because of the gap condition there will be no finite mo- 
mentum divergences either, leading to a finite integral. 
This argument may be extended to models with band 
structure. By "folding out" appropriate bands from the 

first brillouin we may convert the sum X^a / /bz ( ^^ ) 
(where the integral is over the first Brillouin zone) into 



an integral over all of k-space 



//^^M^ 



As 



any possible divergence would come from high energy 
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bands where the dispersion is essentially quadratic and 
the wavefunction is essentially of the continuum model, 



we may reduce the problem to a previously solved case. 
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